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Abstract 



Using a ten dimensional dual string background, we study aspects of the physics 
of finite temperature large N four dimensional SU (N) gauge theory, focusing on the 
dynamics of fundamental quarks in the presence of a background magnetic field. At 
vanishing temperature and magnetic field, the theory has J\f = 2 supersymmetry, 
and the quarks are in hypermultiplet representations. In a previous study, similar 
techniques were used to show that the quark dynamics exhibit spontaneous chiral 
symmetry breaking. In the present work we begin by establishing the non-trivial 
phase structure that results from finite temperature. We observe, for example, that 
above the critical value of the field that generates a chiral condensate spontaneously, 
the meson melting transition disappears, leaving only a discrete spectrum of mesons at 
any temperature. We also compute several thermodynamic properties of the plasma. 



1 Introduction 



In recent years, the understanding of the dynamics of a variety of finite temperature gauge 
theories at strong couphng has been much improved by employing several techniques from 
string theory to capture the physics. The framework is that of holographic [1] gauge/gravity 
duality, in which the physics of a non-trivial ten dimensional string theory background can 
be precisely translated into that of the gauge theory for which the rank (A^) of the gauge 
group is large [21 [3l IH [5], while the number (Nf) of fundamental flavours of quark is small 
compared to N (see ref.[6]). Many aspects of the gauge theory, at strong 't Hooft coupling 
A = QyyiN, become accessible to computation since the string theory background is in a 
regime where the necessary string theory computations are classical or semi-classical, with 
geometries that are weakly curved[2] (characteristic radii in the geometry are set by A). 

These studies are not only of considerable interest in their own right, but have po- 
tential phenomenological applications, since there are reasons to suspect that they are of 
relevance to the dynamics of quark matter in extreme environments such as heavy ion col- 
lision experiments, where the relevant phase seems to be a quark-gluon plasma. While the 
string theory duals of QCD are not known, and will be certainly difficult to obtain compu- 
tational control over (the size of the gauge group there is small, and the number, Nf, 
of quark flavours is comparable to A^) it is expected (and a large and growing literature of 
evidence seems to support this - see below) that there are certain features of the physics 
from these accessible models that may persist to the case of QCD, at least when in a strongly 
coupled plasma phase. Well-studied examples have included various hydrodynamic proper- 
ties, such as the ratio of shear viscosity to entropy, as well as important phenomenological 
properties of the interactions between quarks and quark jets with the plasma. Results from 
these sorts of computations in the string dual language have compared remarkably well with 
QCD phenomenological results from the heavy ion collision experiments at RHIC (see e.g., 
refs. P Uni [m [121 1131 E]); and have proven to be consistent with and supplementary to 
results from the lattice gauge theory approach. 

There are many other phenomena of interest to study in a controllable setting, such as 
confinement, deconfinement (and the transition between them), the spectrum and dynamics 
of baryons and mesons, and spontaneous chiral symmetry breaking. These models provide a 
remarkably clear theoretical laboratory for such physics, as shown for example in some of the 
early work[71 [8] making use of the understanding of the introduction of fundamental quarks. 
Some of the results of these types of studies are also likely to be of interest for studies of 
QCD, while others will help map out the possibilities of what types of physics are available 
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in gauge theories in general, and guide us toward better control of the QCD physics that we 
may be able to probe using gauge/string duals. 

This is the spirit of our current papeiQ Here, we uncover many new results for 
a certain gauge theory at finite temperature and in the presence of a background external 
magnetic field, building on work done recently [T21I2S] on the same theory at zero temperature. 

At vanishing temperature and magnetic field, the large N SU (N) gauge theory has 
M = 2 supersymmetry, and the quarks are in hypermultiplet representations. Nevertheless, 
just as for studies of the even more artificial A/" = 4 pure gauge theory, the physics at finite 
temperature — that of a strongly interacting plasma of quarks and gluons in a variety of 
phases — has a lot to teach us about gauge theory in general, and possibly QCD in particular. 

In section two we describe the holographically dual ten-dimensional geometry and the 
embedding of the probe D7-brane into it. In section three we extract the physics of the probe 
dynamics, using both analytic and numerical techniques. It is there that we deduce the phase 
diagram. In section four we present our computations of various thermodynamic properties 
of the system in various phases, and in sections five and six we present our computations 
and results for the low-lying parts of the spectra of various types of mesons in the theory. 
We conclude with a brief discussion in section seven. 

2 The String Background 

Consider the AdSs-SchwarzschildxS"^ solution given by: 

where dVll = dO"^ + cos^ Odill + sin^ 9d(f)^ , 
and d^ll = d-ip"^ + cos^ tjjdp + sin^ 6*^7^ . 

The dual gauge theory will inherit the time and space coordinates t = and 
respectively. Also, in the solution above, u G [0, 00) is a radial coordinate on the asymp- 
totically AdSs geometry and we are using standard polar coordinates on the S^. The scale 
R determines the gauge theory 't Hooft coupling according to = a' ^J g\y^N . For the 
purpose of our study it will be convenient [Tj to perform the following change of variables: 

=l{u^ + ^^;^^^) = + L\ (2) 
with p = r cos 6 , L = r sin 6 . 
^We note that another group will present results in this area in a paper to appear shortly [25]. 
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The expression for the metric now takes the form: 



( (4r 

Following ref. [6], we introduce fundamental matter into the gauge theory by placing D7- 
brane probes into the dual supergravity background. The probe brane is parametrised by 
the coordinates {xq, Xi, X2, Xs, p, ip, P, 7} with the following ansatz for its embedding: 

(p = const, L = L{p) . 

In order to introduce an external magnetic field, we excite a pure gauge -B-field along the 
directions [T5] : 

B = Hdx^ A (3) 



where if is a real constant. As explained in ref. [T5j, while this does not change the su- 
pergravity background, it has a non-trivial effect on the physics of the probe, which is our 
focus. To study the effects on the probe, let us consider the general (Abelian) DBI action: 



Sdbi = -NjTdi J d'^ dee/\P[Gab + B,,] + 2na'F^f,) , 



(4) 



Ma 



where Tm = Pt/qs = [(27r) ''a'^^f^]"-'^ is the D7-brane tension, P[Gab\ and P[Bah\ are the 
induced metric and induced 5-field on the D7-branes' world-volume. Fab is the world- 
volume gauge field, and Nf = 1 here. It was shown in ref. [U] that, for the AdSs x 
geometry, we can consistently set the gauge field Fat to zero to leading order in a', and 
the same argument applies to the finite temperature case considered here. The resulting 
Lagrangian is: 

,)ll.^'^"^^^±^^tvi.LW. (5) 

16 (p2 + L(p)=) V 1 (61 + 4 (p2 + L(p)2)=) ' ^ " 

For large p^ h, the Lagrangian asymptotes to: 



£^_pV1 + L'(p)2, (6) 
which suggests the following asymptotic behavior for the embedding function L{p): 

L{p) = m + 4 + ■ • • , (7) 
p2 



4 



where the parameters m (the asymptotic separation of the D7 and D3-branes) and c (the 
degree of transverse bending of the D7-brane in the (p, (p) plane) are related to the bare 
quark mass = m/27ra' and the fermionic condensate (ipip) oc — c respectively [H] (this 
calculation is repeated in appendix |A]) . It was shown in ref. [12] that the presence of the 
external magnetic field spontaneously breaks the chiral symmetry of the dual gauge theory 
(it generates a non-zero (ipip) at zero m). However [7], the effect of the finite temperature 
is to melt the mesons and restore the chiral symmetry at zero bare quark mass. Therefore, 
we have two competing processes depending on the magnitudes of the magnetic field H and 
the temperature T = h/ixR^. This suggests an interesting two dimensional phase diagram 
for the system, which we shall study in detail later. 

To proceed, it is convenient to define the following dimensionless parameters: 



P 

m 

This leads to the Lagrangian: 



b 

L{b~p) 



V 



62 



H 



m 



m 



m + ^ + . . . 
p2 



C 



-p' 



1 



1 + 4 p2 + ^(p)2 



1 + L'(p)2. (9) 



For small values of rj, the analysis of the second order, non-linear differential equation for 
L{p) derived from equation ^ follows closely that performed in refs. [71[Tni[IT]. The solu- 
tions split into two classes: the first class are solutions corresponding to embeddings that 
wrap a shrinking in the part of the geometry and (when the vanishes) closes at 
some finite radial distance r above the black hole's horizon, which is located at r = bj\f2. 
These embeddings are referred to as 'Minkowski' embeddings. The second class of solutions 
correspond to embeddings falling into the black hole, since the 5*^ of the Euclidean section, 
on which the D7-branes are wrapped, shrinks away there. These embeddings are referred 
to as 'black hole' embeddings. There is also a critical embedding separating the two classes 
of solutions which has a conical singularity at the horizon, where the 5*^ wrapped by the 
D7-brane shrinks to zero size, along with the . If one calculates the free energy of the 
embeddings, one can show [3, [161 IE] that it is a multi-valued function of the asymptotic 
separation m, which amounts to a first order phase transition of the system (giving a jump 
in the condensate) for some critical bare quark mass mcr. (For fixed mass, we may instead 
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consider this to be a critical temperature.) We show in this paper that the effect of the mag- 
netic field is to decrease this critical mass, and, at some critical magnitude of the parameter 
Tjcv, the critical mass drops to zero. For rj > r^cr the phase transition disappears, and only 
the Minkowski embeddings are stable states in the dual gauge theory, possessing a discrete 
spectrum of states corresponding to quarks and anti-quarks bound into mesons. Further- 
more, at zero bare quark mass, we have a non-zero condensate and the chiral symmetry is 
spontaneously broken. 



3 Properties of the Solution 

3.1 Exact Results at Large Mass 

It is instructive to first study the properties of the solution for m ^ 1. This approxima- 
tion holds for finite temperature, weak magnetic field, and large bare quark mass m, or, 
equivalently, finite bare quark mass m, low temperature, and weak magnetic field. 

In order to analyze the case m 1, let us write L(p) = a + C(p) for a ^ 1 and 
linearize the equation of motion derived from equation (|9]), while leaving only the first two 
leading terms in (p^ + m^)~^. The result is: 



2rf ^ 2{rf + lf-l 



dp{p\') - —f^—.m + zyj_^l_rh + 0(C) = . (10) 



Ignoring the 0{Q terms in equation (10), the general solution takes the form: 

rf . 2(7^2 + 1)2 _i _ 

^^^^ ~ ~ V(m2 + p2)"^ + 96p2(m2 + ^2)3"^ ' (11) 

where we have taken C'(0) = C(0) = 0- By studying the asymptotic behavior of this solution, 
we can extract the following: 

= °-4a5+ 32y +"(y 

= ^ + ^''' + 2''%o(iV (12) 



4a 96a5 \h\ 

By inverting the expression for m, we can express c in terms of fh: 

rf 1 + Arf + 8r]'^ 



4m 96m^ \ m'^ 



+ 0{^] . (13) 
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Finally, after going back to dimensionful parameters, we can see that the theory has devel- 
oped a fermionic condensate: 

iiji;) oc-c = -—H^ + —: . (14) 

Am 9bm-^ 

The results of the above analysis can be trusted only for finite bare quark mass and sufficiently 
low temperature and weak magnetic field. As can be expected, the physically interesting 
properties of the system should be described by the full non-linear equation of motion of 
the D7-brane. To explore these we need to use numerical techniques. 

3.2 Numerical Analysis 

We solve the differential equation derived from equation ^ numerically using Mathematica. 
It is convenient to use infrared initial conditions pTl IT8]. For the Minkowski embeddings, 
based on symmetry arguments, the appropriate initial conditions are: 

L{p)\p=o = L;^, L'(p)|p=o = 0. (15) 
For the black hole embeddings, the following initial conditions: 

L 



Lip)\e.h. = Lin, L'{p)\c.h 

P 



(16) 

e.h. 



ensure regularity of the solution at the event horizon. After solving numerically for L{p) for 
fixed value of the parameter rj, we expand the solution at some numerically large pmax, and, 
using equation ([T]), we generate the plot of — c vs m. It is instructive to begin our analysis 
by revisiting the case with no magnetic field [r] = 0), familiar from refs.[7l [161 lU]- The 
corresponding plot for this case is presented in figure [Tj Also in the figure is a plot of the large 



mass analytic result of equation (13), shown as the thin black curve in the figure, descending 
sharply downwards from above; it can be seen that it is indeed a good approximation for 
rh > fficf Before we proceed with the more general case of non-zero magnetic field, we review 
the techniques employed in ref . [TTj to determine the critical value of fh. In figure |2| we have 
presented the region of the phase transition considerably magnified. Near the critical value 
TTicr, the condensate c is a multi-valued function of m, and we have three competing phases. 
The parameter c is known |8] to be proportional to the first derivative of the free energy of 
the D7-brane, and therefore the area below the curve of the — c vs fh plot is proportional 
to the free energy of the brane. Thus, the phase transition happens where the two shaded 
regions in figure [2] have equal areas; furthermore, for fh < fhcr, the upper-most branch of the 
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Figure 1: The solid curve starting far left (red) represents solutions falling into the black hole, the 
dotted (blue) curve represents solutions with shrinking S^. The vertical dashed line corresponds to 
the critical value of m at which the first order phase transition takes place. The solid black curve 



dropping sharply from above is the function derived in equation (13), corresponding to the large 
mass limit. 



curve corresponds to the stable phase, and the lower-most branch of the curve corresponds 
to a meta-stable phase. For m > rhcr, the lower-most branch of the curve corresponds to the 
stable phase, and the upper-most branch of the curve corresponds to a metastable phase. 
At m = rricr we have a first order phase transition. It should be noted that the intermediate 
branch of the curve corresponds to an unstable phase. 

Now, let us turn on a weak magnetic field. As one can see from figure |3j the effect 
of the magnetic field is to decrease the magnitude of rhcr- In addition, the condensate 
now becomes negative for sufficiently large rh and approaches zero from below as m ^ oo. 
It is also interesting that equation (13) is still a good approximation for m > rricr- For 



sufficiently strong magnetic field, the condensate has only negative values and the critical 
value of rh continues to decrease, as is presented in figure |4j If we further increase the 
magnitude of the magnetic field, some states start having negative values of m, as shown in 
figure [5j The negative values of rh do not mean that we have negative bare quark masses; 
rather, it implies that the D7-brane embeddings have crossed L = at least once. It 
was argued in ref. [7] that such embeddings are not consistent with a holographic gauge 
theory interpretation and are therefore to be considered unphysical. We will adopt this 
interpretation here, therefore taking as physical only the rri > branch of the — c vs rh 
plots. However, the prescription for determining the value of rricr continues to be valid, 
as long as the obtained value of rhcr is positive. Therefore, we will continue to use it in 
order to determine the value of rj = r]„ for which rhcr = 0. As one can see in figure |6l the 
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-c(m) 



0.035 



0.03 



0.025 



m,r:^0.923 




m 



Figure 2: The area below the {—c,m) curve has the interpretation of the free energy of the D7- 
brane; thus the phase transition pattern obeys the "equal area law" — the area of the shaded regions 
is equal. 



-c(m) 



77=0.4 




Figure 3: The effect of the weak magnetic field is to decrease the values of rhcr and the condensate. 



Equation (13) is still a good approximation for m > rhcr- 



value of ?7cr that we obtain is rjcr ~ 7.89. Note also that, for this value of 77, the Minkowski 
m = embedding has a non-zero fermionic condensate €„, and hence the chiral symmetry 
is spontaneously broken. For rj > rj^r, the stable solutions are purely Minkowski embeddings, 
and the first order phase transition disappears; therefore, we have only one class of solutions 
(the blue curve) that exhibit spontaneous chiral symmetry breaking at zero bare quark mass. 
Some black hole embeddings remain meta-stable, but eventually all black hole embeddings 
become unstable for large enough rj. This is confirmed by our study of the meson spectrum, 
which we present in later sections of the paper. The above results can be summarized in 
a single two dimensional phase diagram, which we present in figure [7} The curve separates 
the two phases corresponding to a discrete meson spectrum (light mesons) and a continuous 
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Figure 4: For strong magnetic field the condensate is negative. The value of rricr continues to drop 
as we increase ij. 



-c(m) 




Figure 5: For sufficiently high values of rj there are states with negative m, which are considered 
non-physical. However the equal area law is still valid as long as rricr > 0. 

meson spectrum (melted mesons) respectively. The crossing of the curve is associated with 
the first order phase transition corresponding to the melting of the mesons. If we cross 
the curve along the vertical axis, we have the phase transition described in refs. [3 Uni ITT] . 
Crossing the curve along the horizontal axis corresponds to a transition from unbroken 
to spontaneously broken chiral symmetry |15j. meaning the parameter c jumps from zero to 
Ccr ~ 4.60, resulting in non-zero fermionic condensate of the ground state. It is interesting to 
explore the dependence of the fermionic condensate at zero bare quark mass on the magnetic 
field. From dimensional analysis it follows that: 

cc. = b%,iv) = ~^R'H'/' ■ (17) 

In the T — limit, we should recover the result from ref. P3]: Ccr ~ 0.226R'^H^^^, which 
implies that Ccr{i]) ~ 0.2267]^/^ for 77 ^ 1. The plot of the numerically extracted dependence 
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?7cr^7.89 



-c(in) 




m 



Figure 6: For rj = rjcr the critical parameter rricr vanishes. There are two m = states with equal 
energies, one of them has non-vanishing condensate — Ccr ^ —4.60 and therefore spontaneously 
breaks the chiral symmetry. 

Ccrijf) is presented in figure p| for rj > 77^, Ccxij]) very fast approaches the curve 

0.226773/2. 

This suggests that the value of the chiral symmetry breaking parameter Ccr depends mainly 
on the magnitude of the magnetic field if, and only weakly on the temperature T. 

4 Thermodynamics 

Having understood the phase structure of the system, we now turn to the extraction of 
various of its important thermodynamic quantities. 

4.1 The Free Energy 

Looking at our system from a thermodynamic point of view, we must specify the potential 
characterizing our ensemble. We are fixing the temperature and the magnetic field, and 
hence the appropriate thermodynamic potential density is: 



where (i is the magnetization density and S is the entropy density of the system. Following 
ref . [I9] , we relate the on-shell D7-brane action to the potential density F via: 



dF = —SdT — fidH , 



(18) 



F = 27i^NfTmlD7 , 



(19) 
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m 




leously broken CS 



2 4 6 8 10 12 

Figure 7: The curve separates the two phases corresponding to discrete meson spectrum (light 
mesons) and continuous meson spectrum (melted mesons). 



where (here, Nf = 1): 

pmax 1 



LD7 = U 

Pmiii 



r]=—H; f = r/b; p = p/b; L = L/b; r"^ = + L"^ 



In principle, on the right hand side of equation (19), there should be terms proportional 
to — which subtract the energy of the magnetic field alone; however, as we comment 
below, the regularization of Im is determined up to a boundary term of the form const x H^. 
Therefore, we can omit this term in the definition of F. The boundary action /bound contains 
counterterms designed [20] to cancel the divergent terms coming from the integral in equation 



(20) in the limit of Pmax ^ cxd. A crucial observation is that the finite temperature does not 
introduce new divergences, and we have the usual quartic divergence from the spatial volume 
of the asymptotically AdSs spacetime [21]. The presence of the non-zero external magnetic 
field introduces a new logarithmic divergence, which can be cancelled by introducing the 
following counterterm: 

p4 / n \ r 1 

^1 / f max \ / j4 / -'- 



2 log^^j J d^xV^-B,.B^^ , (21) 

where 7 is the metric of the 4-dimensional surface at p = Pmax- Note that in our case: 

1 



2! 



—lB^,B^^ = H\ (22) 
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20[ 




Figure 8: The solid curves are the numerically extracted dependence Ccr(f?), while the dashed 
curve represents the expected large rj behavior Ccrirj) ~ 0.226r/^/^. The solid curve segments at the 
bottom left and to the upper right (blue) are the stable states. The straight segment and the arc 
that joins it (lower right, red) and red are the unstable states. The rest (cyan) are meta-stable 
states. 

which gives us the freedom to add finite terms of the form const x at no cost to the reg- 
ularized action. This makes the computation of some physical quantities scheme dependent. 



We will discuss this further in subsequent sections. The final form of /bound in equation (20) 
is: 

T 4 d4 i:7'2 1 Pmax /oo\ 

Abound = -^Pmax ~ 2 ' ^ ^ 



It is instructive to evaluate the integral in equation (20) for the L = embedding at zero 
temperature. Going back to dimensionful coordinates we obtain: 

pmax j 

J dppdl + ^ = \pt^y-R^HHogP-^ + ^{l+\ogA-\ogH^ (24) 
' 

The first two terms are removed by the counter terms from /bound, and we are left with: 

F{b = 0,m = 0,H) = 27r'^NfTD7 (1 + log4 - log //2) . (25) 

8 

This result can be used to evaluate the magnetization density of the Yang-Mills plasma at 
zero temperature and zero bare quark mass. Let us proceed by writing down a more general 
expression for the free energy of the system. After adding the regulating terms from /bound, 
we obtain that our free energy is a function of m, b, H: 

F{b, m, H) = 2n^NfTD7b^lD7{rh, rf) + F(0, 0, H) , (26) 
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where iDY^m,!]) is defined via: 
dp \p'(l- 



'1)7 



2^:4 



1 + 



I677 r 



(4f4 + 1)2 



Pmin 

1 



1 



-?7^(1 + log 4 — logr]^ 



V^rn^ - j - plin/4 (27) 



~2 ~2 

r 



In order to verify the consistency of our analysis with our numerical results, we derive an 
analytic expression for the free energy that is valid for m ^ To do this we use that for 



large m, the condensate c is given by equation (13), which we repeat here: 

-2 1 + 4r/2 + 8t]^ 



V 
4m 



96m^ 



+ 0(l/m 



(28) 



as well as the relation dim /dm = —2c. We then have: 

rh 

/I _ 1 + 4^7^ + 81]^ 
c{m,r])dm + ^{rj) = ^(77) - --q^logm ig2^ \-0{l/m^) 



(29) 



where the function ^ (77) can be obtained by evaluating the expression for I^i from equation 
(27) in the approximation L ^ rh. Note that this suggests ignoring the term L'^, which is 
of order c^. Since the leading behavior of at large m is l/m^, this means that the results 
obtained by setting L'^ = can be trusted to the order of 1/m, and therefore we can deduce 
the function ^{rj), corresponding to the zeroth order term. Another observation from earlier 
in this paper is that the leading behavior of the condensate is dominated by the magnetic 
field and therefore we can further simplify equation (27): 



'1)7 





Pmax 


lim 


/ 















^2 

— — log 


- 



(p2 + m^)^ 
2 



- 1 



?7^1ogPmax - -r/M 1 - log 



1 

4 



Comparing to equation (29), we obtain: 

^2 

^iv) = -"§-(3-log^ 
and our final expression for 1 07, valid for m ^ y/fj: 



(30) 



(31) 



'D7 



3 - log 



T] 



1 l±V±^_Lnn/~6^ 
-ri logm hC 1/m 

2 ' ^ 192m4 ^ V / ; 



(32) 



14 



4.2 The Entropy 



Our next goal is to calculate the entropy density of the system. Using our expressions for 
the free energy we can write: 



S 



OF 
df 



H 



ry2dF . ^ ^ ( ~ dim dm dloTdrf 

= ^27r R NfTmb Aim + o— — -- + o—^ — 
GO \ dm do drj'^ do 



-271^ R^NfTorb^ \ Aim + ^cm - 4 



dlpj 2 

Q^2 



2r'R'NfTmb-'S{m,7]' 



(33) 



It is useful to calculate the entropy density at zero bare quark mass and zero fermionic 
condensate. To do this, we need to calculate the free energy density by evaluating the 



integral in equation (27) for L = 0. The expression that we get for Imi^^V^) is: 



1 - 2V1 + -r log 



(34) 



The corresponding expression for the entropy density is: 

1 



S\m=o = 27r''R^NfTmT^ I + a/ 1 + 



WT^ J 



(35) 



One can see that the entropy density is positive and goes to zero as T 0. Our next goal 
is to solve for the entropy density at finite fh for fixed rj. To do so, we have to integrate 



numerically equation (33) and generate a plot of S versus m. However, for m ^ ^ we 



can derive an analytic expression for the entropy. After substituting the expression from 



equation (32) for Im into equation (33) we obtain: 

l + 2rf 



S{m, T] 

or if we go back to dimensionful parameters: 



24m4 



+ 



S{b, m, H) = 27i^R'^NfTmb^ 



b^ + 2R^H^ 
24m4 



(36) 



(37) 



One can see that if we send T 0, while keeping rj fixed we get the T'' behavior described 
in ref. [19], and therefore the (approximate; Nf/N <^ 1) conformal behaviour is restored 
in this limit. In figure jol we present a plot of 5* versus m for r/ = 0.4 . The solid smooth 



black curve corresponds to equation (36). For this S is positive and always a decreasing 



function of m. Hence, the entropy density at fixed bare quark mass m = mb, given by 
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S = 2tt^ R"^ N fT^um^ S / rh^ , is also a decreasing function of m and therefore an increasing 
function of the temperature, except near the phase transition (the previously described 
crossover from black hole to Minkowski embeddings) where an unstable phase appears that 
is characterized by a negative heat capacity. 




m 



Figure 9: A plot of S versus m for ij = 0.4 . The thin (sharply descending and extending to the 



right) black curve corresponds to that large mass result of equation (36). 



4.3 The Magnetization 



Let us consider equation (25) for the free energy density at zero temperature and zero bare 
quark mass. The corresponding magnetization density is given by: 

''°^"0«^'^'''*^'^"?'°^?- '''' 

Note that this result is scheme dependent, because of the freedom to add terms of the form 
const X H'^ to the boundary action that we discussed earlier. However, the value of the 
relative magnetization is given by: 

/i-/io = -(|^) -fio = -2n''R^NfTn7b^l'^] = 2t:^R^N jT^h'^li , (39) 



is scheme independent and is the quantity of interest in the section. In equation (39), 
we have defined jl = —dl£n/dri\m as a dimensionless parameter characterizing the relative 
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magnetization. Details of how the derivative is taken are discussed in appendix [B} The 



expression for yU follows directly from equation (27): 



~^=-^'^ - I ^^~ -y(4,-4 + + le pi + ^^"gP^a. - ^ log -■ (40) 



mm 



For the large m region we use the asymptotic expression for Idj from equation (32) and 
obtain the following analytic result for fl: 

/i = I - I log I + log m + '^^\'^^f^ + 0(l/m^) . (41) 

We evaluate the above integral numerically and generate a plot of jl versus rh. A plot of 
the dimensionless relative magnetization jl versus m for rj = 0.5 is presented in figure [TOj 



The black curve corresponding to equation (41) shows good agreement with the asymptotic 



behavior at large m. It is interesting to verify the equilibrium condition d^/dT > 0. Note 
that since /io does not depend on the temperature, the value of this derivative is a scheme 



independent quantity. From equations (39) and equation (41), one can obtain: 

I - 2^'«^iV,T,,6 (2, - - 2|,) ^ 2.3««iV,T,,^^ > , (42) 

which is valid for large m and weak magnetic field H. Note that the magnetization seems to 
increase with the temperature. Presumably this means that the temperature increases the 
"ionization" of the Yang-Mills plasma of mesons even before the phase transition occurs. 

4.4 The Speed of Sound 

It is interesting to investigate the effect of the magnetic field on the speed of sound in the 
Yang-Mills plasma. Following ref. [19], we use the following definition to thermodynamically 
determine the speed: 

2 S Sd3 + Sd7 , . 

V = — = ■ , 43) 

CV CV3 + CV7 

where cy is the density of the heat capacity at constant volume. To compute the contribution 
coming from the fundamental fiavors in the presence of an external magnetic field, we work 
perturbatively in small Nf/Nc. First, let us recall the adjoint contribution to entropy and 
specific heat [19] : 

Sd3 = ~N^T' , CV3 = 3Sd3 . (44) 
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Figure 10: A plot of the dimensionless relative magnetization jl versus rh for r] = 0.5. The thin 



black curve (starting with a steep descent) corresponds to the large mass result of equation (41). 



To proceed, let us rewrite the entropy density of the fundamental flavours in the following 
form: 



S, 



drf 



(45) 



where M = 2tx'^N fTo-j. The term 4 {Fq — F) /T is simply the contribution from the con- 
formal theory; the deviation from it is related to the conformal symmetry being broken by 
introducing the fundamental flavors and the external magnetic field. This breaking is man- 



ifest by non-vanishing c and r] in equation (45) respectively. Recalling the relation between 



the energy density and the free energy density E = F + TS, and, using equation (45), we 
find that: 



(46) 



Using the definition in equation (43), together with the above results and expanding up to 
first order m. v = Nj/Nc, we get: 



mc m 



l.=.2^c\ XNf7r^7]^(4diD7 , 1 d 



dm J 



(2mc) 



+ 0(z/2). 



(47) 
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The second and third term in equation (47) represent the deviation from the conformal 



value of 1/3 by the presence of the fundamental flavors and the external magnetic field. For 
convenience let us define dv"^ = — 1/3. It is possible to obtain an analytic expression for 
bv^ in the limit of large bare quark mass and small magnetic field (m ^ \f^)- Using our 
previous analytic expressions, we get: 



X- 



N 3 



-T] logm 



-?7^ log 



72m4 



O 



(4J 



It is important to note that equation (48) is valid only up to first order in z/. To proceed 



beyond the large bare quark mass and small magnetic field limit, we study numerically the 



velocity deviation, which is summarised in figure 11, We observe from figure [11(a) that, for 
small magnetic field, the deviation is similar to the zero magnetic field case; 6v'^ approaches 
zero (corresponding to restoration of the conformal symmetry) from below in both the T — 
and T — > cxD limits. However in presence of large magnetic field (see figure 11(d)), we see 
that > 0, and the conformal value is never attained. 



5 Meson Spectrum 

In this section, we calculate the meson spectrum of the gauge theory. The mesons we are 
considering are formed from quark-antiquark pairs, so the relevant objects to consider are 
7-7 strings. In our supergravity description, these strings are described by fluctuations (to 
second order in a') of the probe branes' action about the classical embeddings we found 
in the previous sections[22j. Studying the meson spectrum serves two purposes. First, 
tachyons in the meson spectrum from fluctuations of the classical embeddings indicate the 
instability of the embedding. Second, a massless meson satisfying a Gell-Mann-Oakes-Renner 
(GMOR) relation will confirm that spontaneous chiral symmetry breaking has occurred. As 
a reminder, in ref. [22], the exact meson spectrum for the AdSs x background was found 
to be given by: 

2777 / 

M{nJ) = —^{n + e + l){n + i + 2) , (49) 
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Figure 11: The deviation of the speed of sound from the conformal value in units of (z^A)7r^/3 in 
the presence of fundamental matter and an external magnetic field. The curves coming in from 
the left (red) correspond to black hole embeddings, and the curves coming in from the right (blue) 
correspond to Minkowski embeddings. The vertical dashed line represents the phase transition 
point, and the flatter dashed curves (green) correspond to the approximate analytic expression 



given in (48). We do not include the curve of the analytic result in 11(d) since the approximate 
formula is not valid for high magnetic fields. 



where £ labels the order of the spherical harmonic expansion, and n is a positive integer that 
represents the order of the mode. The relevant pieces of the action to second order in a' are: 



S/Nf 



(4) 



-Td7 j d^i^Qab + Bab + ^Tia'Fab + (2^') ^l^ 

+ {2Tia' f firl [ F(2) A F(2) A P [^(4)] , 



Ms 



F(2) A 5(2) A P 



(4) 



1 



dt A dx^ A dx^ A dx'^ , 



(7(4) = (1 — cos'* 6*) sin i}) cos dtp A (i02 A dcp^ A 

9s 



(50) 
(51) 

(52) 
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where P [C(4)] is the pull-back of the 4-form potential sourced by the stack of A^^^ D3-branes, 
P C*(4) is the pull-back of the 4-form magnetic dual to C(4), and -F(2) is the Maxwell 2- 
form on the D7-brane worldvolume. At this point, we resort to a different set of coordinates 
than we have been using. Instead of using the coordinates (p, L) introduced in equation ([2]), 
we return to the coordinates [z = l/u'^,6) because the analysis is simpler. We consider 
fluctuations of the form: 



9 = 9o{z) + 2na'xin 

h = 2w<i>(r) , 



(53) 
(54) 



where the indices a, 6 = ... 7 run along the worldvolume of the D7-brane. 9q{z) corresponds 
to the classical embedding from the classical equations of motion. Plugging the ansatz in 
equations (53) and (54) into the action and expanding to second order in (27ra'), we get as 
second order terms in the lagrangian: 
1 



— i2<J)2 

— Cp2 

/-WZ 



+ 2^ 



9|v^ - E'^RXdeV^ 



V^S^^R^ sin^ 9oda^db^ 



4 

XF23 



XF23f 



-^FoiBisR^ sin ^ COS tfjd, (l - cos^^o) = -$^01^23^^ sin ^ cos V'^^i^ -(55) 



We have taken Eat = g^^b + ^ab to be the zeroth order contribution from the DBI action. 
In addition, we use that E"^ = S"^ + J''^ where S"'' = S'''' and J"'' = -J^". We use this 
notation for brevity. The indices m,n,o,p = 4 ... 9 run in the transverse directions to the 
D3-branes. From these lagrangian terms, we derive the equation of motion for x to be: 



= dj ^^S^^R^ 



Ah'^z 



,4^4 



23 



1 + 4^2 ( 

d, (y/^R'^9'QE''J^^ 



The equation of motion for $ is given by: 



da [y/-ES''^R^ sin^ ^0^6$) - ^^01^23^^ sin cos tl)d,K = . 



(56) 



(57) 
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The equation of motion for Aj, is given by: 

da (^-^S'^'S'^'Fa^,, - Xf (^2^3 ' ^S^^ + B.s^d^K {S^.S^ - S^) 

I 1 ^ mnopsra sjb p 



= . 



(58) 



We are allowed to set A„i = with the constraint (using that 5*^^ = 5*^^): 

Therefore, we can consistently take Aq — diAi — 0, ^2^2 = —d^A^. With this particular 
choice, we have as equations of motion for the gauge field: 

-do (^^/^S'^'S^'doAi^ + d,KB23do^ - (^V^S''S''d,A,''^ 

-drniy^S'^'^S^'dnAi] = 0, 



-do (x/=:B5°'^522aoA2) + fd,x - [^s/^S^'S^^d^A^ 

-drn (y^S'^^'S^^dnA^ 

-do {^s/^S'^'S'^doAs^ - fd2X - (^^/^S^'S^'d^As 

-d^ ( V^S'^^'S^'dnAs 



0, 



= 0, 

where the indices rh, n run over the that the D7-brane wraps. If we assume that diX — 0, 
we find that the equations for A2 and A3 decouple from x- Therefore, we can consistently 
take F23 = 0, or, in other words, A2 — A3 — 0. This simplifies the equations of motion that 
we need to consider to: 



= da ^s"^rU ^^"'^^ \d,x -X Ql^f^ 

-d,{E''R^d'od6^r^^ , 
= -do [y^S'^^'S^^doA^ + d,KB2zdo^ - d, [sf^S'^S^^d.A^ 
-d^iy^S'^'^S^^dr^A^ , 



(59) 



(60) 

= da{y^S''^Pi^^m^eodh^-FoxB2zR^?^\^i)co^i)d^K . (61) 

In the proceeding sections, we will work out the solutions to these equations numericahy 
using a shooting method. With an appropriate choice of initial conditions at the event 
horizon, which we explain below, we numerically solve these equations as an initial condition 
problem in Mathematica. Therefore, the D.E. solver routine "shoots" towards the boundary 
of the problem, and we extract the necessary data at the boundary. 



22 



5.1 The X Meson Spectrum 



In order to solve for the meson spectrum given by equation (59), we consider an ansatz for 
the field x of the form: 

X = h{z) exp {—iut) , (62) 

where we are using the same dimensionless coordinates as before, with the addition that: 

z = b~'^z , u = Cj . 

In these coordinates, the event horizon is located aX z = 1. Since there are two different 
types of embeddings, we analyze each case separately. We begin by considering black hole 
embeddings. In order to find the appropriate infrared initial conditions for the shooting 
method we use, we would like to understand the behavior of h{z) near the horizon. The 
equation of motion in the limit of 5 ^ 1 reduces to: 

h"{z) + -J—h'{2) + — Mz) = . (63) 

-2- 1 16(5- 1)^ 

The equation has solutions of the form (1 — z)^*'^/^, exactly of the form of quasinormal modes 
[23] . Since the appropriate fluctuation modes are in-falling modes [21], we require only the 
solution of the form (1 — z)~*'^/^. This is our initial condition at the event horizon for our 
shooting method. In order to achieve this, we redefine our fields as follows: 

h{z) = yCz)il-~zr^/\ 

which then provides us with the following initial condition: 

y{z~^l) = e, (64) 

where e is chosen to be vanishingly small in our numerical analysis. The boundary condition 
on y'{z 1) is determined from requiring the equation of motion to be regular at the 
event horizon. The solution for the fluctuation field y{z) must be comprised of only a 
normalizable mode, which in turn determines the correct value for uj. Since we are dealing 
with quasinormal modes for the black hole embeddings, uj will be complex; the real part 
of UJ corresponds to the mass of the meson before it melts, and the imaginary part of uj is 
the inverse lifetime (to a factor of 2) [24j. We begin by considering the trivial embedding 
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6q{z) = (a black hole embedding). This embedding corresponds to having a zero bare 



quark mass. The equation of motion (59|) simplifies tremendously in this case: 

22 1 



h"{2) + 



P-1 2(1 



z'^rf) 



42:2 (^2 _ lY 



. 



(65) 



We show solutions for u in figure 12 as a function of the magnetic field r/. In particular, 
we find the same additional mode discussed in ref. [25]. This mode becomes massless and 
eventually tachyonic at approximately r] ^ 9.24. This point was originally presented in 
figure [8} where the intermediate unstable phase joins the trivial embedding. This is exactly 
when the — c vs fh plot has negative slope for all black hole embeddings. 
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Figure 12: The x meson mass as a function of magnetic field for the trivial embedding. The 
upper (red) curve is the generahzation of the mode described in ref. |24j. The lower (blue) curve is 
the generalization the mode discussed in ref. [25 . We do not extend this second curve to small r] 
because the numerics become unreliable. 

We now consider embeddings with non-zero bare quark mass. This means solving the full 



equation (59). We have both embeddings to consider; for the black hole embeddings, we 
will follow the same procedure presented above to solve for the complex oj. We can still use 
the same procedure because in the limit of z -h^ 1, the equation of motion still reduces to 



equation (63). For the Minkowski embeddings, we do not have quasinormal modes, and Cj 



is purely real. Therefore, we use as initial conditions: 



X{z Zmax) = e , (66) 

X{Z Zmax) = OO . (67) 
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In figures 13 and 14, we show solutions for very different magnetic field values. In the former 
case, rj is small, and we do not have chiral symmetry breaking; in the latter, t] is large, and 
we have chiral symmetry breaking. It is important to note that in neither of the graphs do 



we find a massless mode at zero bare quark mass. In figure 13, we find that fluctuations 
about both the black hole and Minkowski embeddings become massless and tachyonic (we 
do not show this in the graph). The tachyonic phase corresponds exactly to the regions in 
the — c vs m plot with negative slope. 




Figure 13: The x meson mass as a function of bare quark mass for r] = 1. The dashed (blue) 
curve corresponds to fluctuations about black hole embeddings. The solid (red) line corresponds to 
fluctuations about Minkowski emeddings. These modes have a purely real uj. The straight dashed 
(black) line corresponds to the pure AdSs x solution. 



6 The <l> and A Meson Spectra 



Let us now consider the coupled fluctuations of $ and A in equations (60) and (61). We 
consider an ansatz (as before) of the form: 



$ = (j){z) exp {—iuot) 
Ai = A{z) exp {—iut) 
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Re[w] = a> 




Figure 14: The x meson mass as a function of bare quark mass for r] = 10. The dashed (black) 
line corresponds to the pure AdSs x solution. 



It is interesting to note that, for the trivial embedding 6q{z) = 0, one of the coupled equations 
is equal to zero, and we simply have: 



^ ' [z^ - 1) {I + z^T]^) ^ ' Miz'^-Xf ^ ' 



Again, we note that in the limit of 2; — >^ 1, we have: 



A!\~z) + ^A{z)+ '^'^ = 

z-\ 16(5-1)^ 



This is exactly the form of equation (63), so A{z) has the same solutions of in-falling and 
outgoing solutions, which provides us with the necessary initial conditions for our shooting 
method. We show solutions for a) in figure [TSj Unfortunately, we do not know how to solve 
for the quasinormal modes for other black hole embeddings. Since the equations of motion 
are coupled, we are unable to find an analytic solution for the fluctuations near the event 
horizon. This prevents us from using infrared initial conditions for our shooting method. 
However, we are actually more interested in searching for the "pion" of our system, which 
will occur when we have chiral symmetry breaking. In those cases, we are only dealing with 
Minkowski embeddings with a pure real cj, so we may ignore the black hole embeddings. 
Since the equations of motion are coupled, it turns out that only for specific initial conditions 
will both fluctuations only be comprised of normalizable modes. We represent this by a 
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Re[w] 




Figure 15: The A meson mass as a function of magnetic field for the trivial embedding. 

parameter a (which we must tune) as follows: 

A(5 ^ 5max) = icosa , (68) 
0(5 -> 5max) = sin a . (69) 

The initial conditions on A' and 0' are determined from the equations of motion. We show 



several solutions in figure 16 There are several important points to notice. First, we find 



that the lowest mode satisfies an GMOR relationship given by: 

« l.lm^/2 . (70) 

Therefore, we have found the Goldstone boson of our system related to the breaking of chiral 
symmetry. Second, the modes exhibit the Zeemann splitting behavior that was discussed in 
ref. [15]. It is interesting that the various modes always cross each other at approximately 
m ^ 2. We have no intuitive explanation for this behavior. 



7 Conclusions 

We have extended the holographic study of large gauge theory in an external magnetic 
field, started in ref. [15] , to the case of finite temperature, allowing us to study the properties 
of the quark dynamics when the theory is in the deconfined plasma phase. 

The meson melting phase transition exists only below a critical value of the applied 
field. This is the critical value above which spontaneous chiral symmetry breaking is triggered 
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Figure 16: Coupled A — (j) fluctuations for rj = 10. Tlie daslied (blaclt) line corresponds to the pure 
AdS5 X solution. 

(in the case of zero mass). Above this value, regardless of the quark mass (or for fixed quark 
mass, regardless of the temperature) the system remains in a phase with a discrete spectrum 
of stable masses. Evidently, for these values of the field, it is magnetically favourable for the 
quarks and anti-quarks to bind together, reducing the degrees of freedom of the system , as 
can be seen from our computation of the entropy. Meanwhile, the magnetization and speed 
of sound are greater in this un-melted phase. 

There have been non-perturbative studies of fermionic models in background mag- 
netic field before, and there is a large literature (see e.g., the reviews of refs. [271 [2H], and 
the discussion of ref.[29] and references therein). Generally, those works use quite different 
methods to examine aspects of the physics — some primary non-perturbative tools are the 
Dyson-Schwinger equations in various truncations). Our results (and the zero temperature 
result obtained with these methods in the zero temperature case [15]) are consistent with the 
general expectations from those works, which is that strong magnetic fields are generically 
expected to be a catalyst for spontaneous chiral symmetry breaking in a wide class of models 
(see e.g., refs. [3U1 [221 12Z] for a discussion of the conjectured universality of this result). 

While it is satisfying that our supergravity/string methods, which probe the gauge 
theory non-perturbatively via the holographic duality, confirm those other non-perturbative 
approaches, it would be interesting and potentially useful to compare the results in more 
detail, as this would (for example) allow a better understanding of the systematics of the 
Dyson-Schwinger truncation schemes. Whether or not such a direct comparison of these 
very different non-perturbative approaches is possible would also be interesting to study 
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in its own right, potentially shedding light on other non-perturbative phenomena in field 
theory that have been studied using such methods. This avenue of investigation is beyond 
the scope of this paper, however, and we leave it for future study. 
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A Calculating the Physical Condensate 

The condensate (density) is given by: 



SF 



6mn 



(71) 



where F is the free energy density and ruq is the physical bare quark mass. The free energy 
density is given by equation (26), which we duplicate here: 

1 



2n'NfTD7b 



pmax 



Lpn 



1 



r;4 



1 



2^4 



1 



IGf] r 



(4f4 + ly 



1 + L'2 - p' 



;;3 



^ Pmin 2 ■ 



-?7MogPmax - o?7^(l + log 4 - log?7^ 
) 

Therefore, to calculate the condensate, we have: 



6F ,5F 2-na' 5F 2W 5F 

2-na 



5mq Sm b Sm b 5L(e) 



(72) 



where we are using a new set of coordinates z = 1/p such that e = 1/pmax- In order to 
continue, we consider the variation of the free energy density: 



6F = 27i^NfTD7b'^ 



dz 



dC 



d dc 



dL (5) dz dV (z) 



6L (z) + -^5L ii) 



.(73) 



The first term is set to zero by the equation of motion. The boundary term evaluated at 
^max is zero since for Minkowski embeddings L'(5max) = and for black hole embeddings 
1 — l/16f® = 0. Therefore, we are left with: 



5F 



-2n'NfTn7b'-y—SLie) 
oL' (z) 



-ATT^NfTDrb^cSTJi 



(74) 



where we have used the asymptotic expansion of L(e) 
expression for the condensate is given by: 



m 



+ cz"^. Therefore, our final 



-S-K^a'NfTmb^c = , c 



(75) 
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B The 77 Variation 



In calculating the magnetization, we need to consider the on-shell quantity: 



61] 



; Id7= I dpC{p,r];L,L') . (76) 

T 



In our discussion, fixing the temperature is equivalent to fixing the bare quark mass. This 
is true because the two quantities (T, m^) are related inversely to each other by the dimen- 
sionless quantity rh = m/b = 2TTa'mq/TTR^T . Therefore, we are interested in calculating: 

dpi---— + -~ — + ^] = [---—] + dp— , (77) 



Sf] I J yOL' 6r] dL Sr] dr] I \ dU St] I _ J dr] 
where we have used the equation of motion to simplify the expression. We can further 



simplify the expression by noting that the boundary term in equation (77) is zero, because 



=0, and — L = — = 0. 78 

The last relation follows from the fact that we are taking the variational derivative with 
respect to 77 at fixed in. Therefore, the necessary computation is simplified tremendously to: 

' ^ \ / ~ \ Pmax ^ 

/ T \ / m . 
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